Discrete scale invariance and stochastic Loewner evolution 
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In complex systems with fractal properties the scale invariance has an important rule to classify 
different statistical properties. In two dimensions the Loewner equation can classify all the fractal 
curves. Using the Weierstrass-Mandelbrot (WM) function as the drift of the Loewner equation we 
introduce a large class of fractal curves with discrete scale invariance (DSI). We show that the fractal 
dimension of the curves can be extracted from the diffusion coefficient of the trend of the variance 
of the WM function. We argue that, up to the fractal dimension calculations, all the WM functions 
follow the behavior of the corresponding Brownian motion. Our study opens a way to classify all 
the fractal curves with DSI. In particular, we investigate the contour lines of 2D WM function as a 
physical candidate for our new stochastic curves. 



INTRODUCTION 



Discrete scale invariance (DSI) is a weaker kind of scale 
invariance according to which the system obeys scale in- 
variance 0{\x) = pi(X)0(x) only for specific choices of A. 
The list of applications of discrete scale invariance (DSI) 
covers phenomena in programming and number theory, 
diffusion in anisotropic quenched random lattices, growth 
processes and rupture, quenched disordered systems, tur- 
bulence, cosmic lacunarity, etc, for review of the concept 
and applications see 0, 0| . In some cases the DSI is man- 
ifested in the geometry of the system, the most famous 
cases being animals Q and Diffusion-limited aggregation 
(DLA) [4] . The fractal dimension of these models shows 
log-periodic corrections to scaling which is the signature 
of the presence of complex exponents and DSI 
One of the methods to tackle these models, especially 
the DLA, is the iterative conformal mapping technique 
introduced by Hastings and Levitov in |7|. It is based 
on approximating an addition of a particle to preexist- 
ing domain with bumps at nearby positions produced by 
conformal maps. This method was used successfully to 
describe the Laplacian stochastic growth models [7| as 
well as the non-Laplacian transport processes [$], for a 
recent review see [9(. 

Another method, in the same spirit as the Hastings and 
Levitov method, to study geometrical critical systems is 
by Schramm-Loewner evolution 



[lOj . for review see [llj. 
It is a family of random planar curves with conformal 
symmetry that can be described by successive random 
conformal maps. This method recently has found lots of 
applications in the physical systems such as domain walls 
in statistical mode ls [lllj . zero-vorticity lines of Navier- 
Stokes turbulence [12| and iso-height lines in rough sur- 



faces [Ij,ll4|. 



Although SLE gives a powerful way to classify all the 
conformally invariant curves in two dimensions it is still 
too restrictive to classify the curves with weaker sym- 
metries such as scale invariance. In this letter we will 
introduce a broader range of fractal curves in two dimen- 
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FIG. 1: (Color on line) Curves, from left to write k = 
0.25, 1,2.5, generated by taking WM function with A = 2.5. 
The total time is T = 10 and the number of points is 
N — 3 x 10 4 . The computer running time for every real- 
ization is 15 min. 



sions. The method takes the advantage of the Brownian 
motion limit of the WM function and opens a way to 
classify the DSI fractal curves in two dimensions. 
The idea behind SLE formulation is parameterizing a 
growing curve with time t in a two dimensional domain 
with some successive conformal maps that can remove the 
curve. The most familiar case is the curve growing with 
the tip at j t in the upper half plane H. The curve can be 
parametrized with the conformal map gt(z) : M \ Ift — > H 
which maps the upper half plane minus the hull [15| K t , 
to the upper half plane and obeys the following Loewner 
equation 



dg t {z) 



2dt 



(1) 



g t (z)-U(ty 

where Ut is the drift of the process and can be an arbi 
trary continuous function with Holder exponent h > , 
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Every growing curve is related to a special U(t) and by 



can conclude that gt{z) satisfies the same Loewner equa- 
tion as ([T]). Using this property one can claim that gt(z) 
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FIG. 2: (Color on line) The fractal dimension of the gen- 
erated curves versus re for different drifts: a) ordinary SLE 
with Brownian motion as a drift, b) W SLE(e tx , 2.5, re) and c) 
WSLE(sin~ 1 (sin(x)), 2.5, re). The continuous lines are the 
lines coming from the equation (J6|. In a the number of real- 
izations for every re is 1000 and the error bares are smaller 
than the diameter of the circles. In b and c the number of re- 
alizations for every re is 500. The number of points on every 
realization was N = 3 x 10 4 . 

having U(t) one can extract a growing curve in the upper 
half plane. Considering Ut as a process proportional to 
the Brownian motion ^JTiBt one can extract conformal in- 
variant curves and the evolution is called SLE ■ Brow- 
nian motion is a Gaussian Markov stochastic process with 
mean zero and E[BtB s ] — min(t,s). All conformal in- 
variant curves in two dimensions can be parametrized by 
a Brownian motion. 

The fractal dimensions of the curves are linearly propor- 
tional to re and given by 
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FIG. 3: (Color on line) The fractal dimension of the generated 
curves versus A for two different res for W SLE(e" , A, l)(top) 
and WSLE(e lx , A, 0.5) (bottom). The number of realizations 
for every A is 50. Inset: The Lomb periodogram of the Re(7t), 
here u) n is equal to A n . 

is the conformal map that takes the complement of the 
hull K t :— K\2 t /X and maps it to the upper half plane. 
The natural consequence of taking discrete scale invariant 
drift in the Loewner equation is getting a discrete scale 
invariant curve in the upper half plane and vise versa. 
With this introduction, in this paper we will take the 
generic process with DSI, Weierstrass-Mandelbrot (WM) 
type functions 17 -Til, as a drift and we will investigate 
numerically the fractal properties of the corresponding 
curves in the upper half plane. We will also introduce 
the contour lines of the 2D WM functions as the new 
fractal curves with DSI and as the possible candidates 
for our new stochastic Loewner evolution. 



(2) 



The linearity of the relation between fractal dimension 
and re comes from the specific form of the Fokkcr-Planck 
equation of the process which is a second order differen- 
tial equation. The natural question is what is the ap- 
propriate U(t) for the curves with DSI? To answer this 
question consider a drift with DSI as U{t) := jU(X 2 t), 
where A is a positive real number, then one can show that 
g t (z) = x9xh(^ z ) satisfies the same Loewner equation as 
(UJ). To see this consider the new time parametrization in 
equation (fT]) as t —¥ a(t). Then one can write the equa- 

2 da(t) 

tion yj) for the general drift as dg t (z) — jjTj^r^m ; 
where gt{z) = g a (t){ z ) an d go(z) = z. Considering 
a(t) = X 2 t it is easy to get 



LOEWNER EQUATION WITH DSI DRIFT 



dg t {z) 



where gt(z) = and go{z) = z. Since U(t) is the same 
as U(i), in the distribution sense, for our DSI process one 



2dt 



g t (z)-{U(XHY 



(3) 



The building block of the functions with DSI is 
Weierstrass-Mandelbrot (WM) function which was stud- 
ied by Berry and Lewis in [171 ■ It i s a random contin- 
uous non diffcrentiable mono fractal function that one 
can see it in a wide range of physical phenomena such 
as sediment and turbulence [l|, fractal properties of 
the landscapes and other enviromental data [201] , contact 
analysis of elastic-plastic fractal surfaces [21( and prop- 
agation and localization of waves in fractal media [221 ] . 
The WM function can be defined as follows 



W(h,t)= J2 (HO) - h(X n t))e^ 



(4) 



where h(t) can be any periodic function which is differ- 
entiable at zero and <p n is an arbitrary phase chosen from 
the interval [0, 2tt], to make the series convergent we need 
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to consider < H < 1. It was proved in 17|, [23J that the 
above process with h(x) = e 11 converges to the fractional 
Brownian motion in the limit A — > 1 . In our numerical 
calculation we will take h(x) = e lx , sin _1 (sin(a;)), how- 
ever, our results are extendable to the most generic cases. 
The WM process is not Markov but has stationary incre- 
ments. It is a fractal function for the finite n m i n and 
also for n m i n — — oo, however, it has full DSI just for 
n m in — — oo. The fractal dimension of the function is 
2 — H where H is the Holder exponent. 

Taking Re(y/KW(h, t)) as the drift of the Loewner 
equation one can get simple curves just for H > h. The 
fractal dimension of the curves is 1 for H > | but non- 
trivial for H = 7}. From this moment, we will consider 
just H = i which has the same fractal dimension as the 
Brownian motion but it is not a Markov process. One 
can extract the corresponding curves numerically as an 
iteration process of infinitesimal conformal mappings as 
follows: y(j't) = fi°f2 ---°fj(£,j), where £j is the drift, in 
our case WM process, and fj(z) 



is the inverse Loewner map that can produce a small slit 
at £ 3 . The discretized driving force is constant in the in- 
terval [(j — 1)t, jr] where r = jj is the mesh of the time 
and N is the number of points on the curve. We will 
call these new curves WSLE(h, A, k). In Fig 1 one can 
see WSLE(e lx , 2.5, n) for three different k's. In Fig 2 
the box fractal dimension of SLE, WSLE(e lx , 2.5, n) and 
W S LE (sin" 1 (sin(x)) , 2.5, k) were shown with respect to 
k and we observe that they are all linearly dependent on 
k. The slopes of the curves in Fig 2 is dependent on the 
function h(x) and A. 

One can find a theoretical approximation for the fractal 
dimension of the extracted curves by calculating the sec- 
ond moment of the real part of the WM function. Using 
Poisson re-summation formula one can find the following 
formula for the trend of the WM process 



\Re(^W(h,t + T)) -Re(^W{h,t))\ 2 
kt r°° h(0)-h(x) ]2 
lo x 



2 In A 



dx. 



(5) 



The above equation shows that Vln A Ke(W(h, t)) con- 
verges to the Brownian motion in the limit A — > 1. Then 
if we consider the WM process as a non-Markovian per- 
turbation of the Brownian motion, one can introduce an 
effective diffusion constant 



1 



2 In A 



h(0)-h{x) 



'dx, 



(0) 



Assuming the universality of the coefficient g the fractal 
dimension will be the same as ([2]) but with the k sub- 
stituted with K e ff. It turns out that the above approxi- 
mation gives convincing results, in Fig 2 the continuous 
lines are the lines coming from this approximation. To 
check the relation between fractal dimension and A we 
calculated the fractal dimension for k = 0.5 and 1 for 




FIG. 4: (Color on line) The contour lines of the 2D WM 
function with A = 1.5, M = 20 and L = 200. 



the different As, the results are shown in Fig 3. Up to 
the numerical accuracy the proposed formula is perfectly 
confirmed. This result shows that as long as the fractal 
dimension calculations is concerned, the Brownian mo- 
tion behaves like a fixed point of the all the DSI processes 
with H = i. In the same figure we put also the result 



of the Lomb test [24( on the Re(7t) which shows the DSI 
properties of the curves. 

In the rest of the paper we will introduce some new 
DSI curves that can be studied with WSLE(h, A, k). 



CONTOUR LINES OF THE DSI SURFACES 

It was shown in [l4| that the contour lines of Gaussian 
free field theory (GFF) can be described by SLE(n,p). 
In other words the contour lines of rough surfaces with 
the zero roughness exponent are conformally invariant 
and can be described by the Schramm's method. Having 
this fact as a motivation we generated a 2D WM function 

m 

In A ^ °° 2ir\ n 

m— 1 n— — oo 



.(22^cos(fl-a w ) + ^ n , n )),(7) 



1/2 

where p = (a; 2 + j/ 2 ) and 8 = tan (^). The phase 
a m = nm/M and A m make the system homogeneous 
and ^ mjrl , which is randomly chosen from the interval 
[0, 27r], makes the system random. For the isotropic sur- 
faces we have A m = const. In numerical methods the 
two parameters M and n and the size of the system L 
should be finite but large enough that z (x, y) does not 
change with increasing these parameters. The roughness 
exponent and the fractal dimension D = 3 — H are two 
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FIG. 5: (Color on line) The perimeter of the contour lines 
I with respect to the radius of gyration R. To get the graph 
we averaged over all the perimeters of the contour lines with 
the radius of gyration in the interval [ln(i?), ln(i?) +<5]. Inset: 
The Lomb periodogram of the contour lines. 



important parameters in simulating WM function. The 
surfaces with larger H seems smoother than the surfaces 
with smaller H . The 2D WM function obeys the scaling 
law z(Xx, Ay) = X H z(x,y). This shows that z(x,y) is a 
self similar function with respect to the scale A. 

We generated 750 realizations of the above surfaces 
with H = 0, L = 1024, M = 20 and A = 1.5 then we 
constructed the contour lines of the surfaces (as already- 



explained in [2Q, |27|) at the average height, see Fig. 4. 



The average fractal dimension of the roughly 10 6 contour 
lines is D — 1.505 ± 0.005, see Fig 5, which is very close 
to |. In the next step we took the largest contour lines 
that do not cross the boundary and found |fj|, where 
= — R| is the distance of the points of the contour 
r, from the center of mass R. This parameter shows very 
clear DSI after using Lomb periodogram, Fig 5. 
To check the DSI of the drift of the Loewner equation 
directly, we first consider an arbitrary placed horizon- 
tal line representing the real axis in the complex plane 
across the 2D WM profiles. Then we cut the portion of 
each curve "ft above the real line as it is in the upper 
half plane H. There are different numerical algorithms 
to extract the driving function of a given random curve 
by inversion of the Loewner equation [28| . We have used 
the successive discrete, conformal slit maps based on the 
piecewise constant approximation of the driving function 
that swallow one segment of the curve at each time step 
29}. To do this, at first, we have determined the con- 
tour lines of the 2D WM surfaces as described above 
by the sequences of points Zq, z%, Zjq. We obtained 
2000 such curves with an average number of points about 
5000. After that all of the contour lines were mapped by 
tp(z) = znz/(zn — z). To avoid numerical errors only the 
part of the curves corresponding to capacity 8 — 300 were 
used. It is worth to mention here that since we do not 
expect to have conformal invariance for our contour lines 
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FIG. 6: (Color on line) The variance of the drift of the 
Loewner equation for the contour lines of 2D WM function. 
Top inset: the Lomb periodogram of the drift. Bottom inset: 
the correlation of the increments of the drift of the contour 
lines of the GFF and the contour lines of the 2D WM function. 
The number of realizations were approximately 7000 and the 
number of points on every realization was 3000 — 4000. 



the effect of the map f{z) on the measure of the curves 
is not clear. Since this map is a global conformal map 
we do not expect a huge difference between the measure 
of curves before and after the mapping. Fig. 6 shows the 
properties of the drift calculated with the above method. 
The drift shows clear DSI in the Lomb periodogram with 
A = 1.47 ± 0.09, moreover the variance of the drift is lin- 
ear around the line with slope K — 3.9 ± 0.2 which, up to 
the numerical errors, is very close to 4, see Fig 6. The 
increments of the extracted drift shows oscillating non- 
zero correlation which is clearly indicating that the drift 
is not Markovian, in Fig 6 we compared these results 
with those coming from the simulation of the contour 
lines of the Gaussian free field theory. This comparison 
shows that the contour lines of 2D Weierstrass function 
with H = could be related to the Loewner equation 
with the WM function as a drift. Here it is worth men- 
tioning that we expect the same results for the 2D WM 
functions generated by the other periodic functions in- 
steed of cos(x), we will discuss properties of the contour 
lines of 2D WM function for generic H for the different 
periodic functions in the forthcoming paper [3p| ] . 



CONCLUSION 

We made the first attempt to classify and study sys- 
tematically the fractal curves with DSI based on Loewner 
evolution. Our approach is from many sense reminiscent 
of the renormalization group and universality ideas in 
statistical mechanics and field theory. By looking to the 
Brownian motion as the fixed point of all the WM func- 
tions with H — \ we showed that the fractal dimension 
of the fractal curves follows the behavior of the fixed 
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point. Our approach supports the idea of looking to the 
scale invariant curves as the perturbation of the confor- 
mally invariant curves. We also extracted contour lines 
of 2D fractal surfaces as physical examples for the fractal 
curves with DSL In addition the corresponding drift of 
the Loewner equation for these curves was extracted by 
using numerical methods. The drift shares many proper- 
ties with the WM process. 

It is interesting to understand the connection of our ap- 
proach to the older studies bases on iterative conformal 
mappings [7| , see also [31| . The DSI is manifested mostly 
because of the hierarchical underlying graphs, our study 
shows that it may be possible to study statistical models, 
on graphs with the well-defined continuum limit, which 
show DSI at the critical point. Investigating statistical 
models with such a property could be interesting. 
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